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ABSTRACT 

We present a new zoom-in hydrodynamical simulation, “Eris5i7 ” which follows the 
cosmological evolution and feedback effects of a supermassive black hole at the center of 
a Milky Way-type galaxy. EiisBH shares the same initial conditions, resolution, recipes 
of gas cooling, star formation and feedback, as the close Milky Way-analog “Eris”, but it 
also includes prescriptions for the formation, growth and feedback of supermassive black 
holes. The aim of this simulation is to study the evolution of supermassive black holes 
and AGN feedback effects in a late-type galaxy. We find that the galaxy’s central black 
hole grows mainly through mergers with other black holes coming from infalling satellite 
galaxies. The growth by gas accretion is minimal because very little gas reaches the sub- 
kiloparsec scales, possibly a reflection of the fact that spiral galaxies have a long period of 
quiescent evolution in which their potential is only weakly perturbed by external triggers. 
The final black hole is, at z = 0, about 2.6 million solar masses and it sits closely to the 
position of SgrA* on the Mbh — ^Buige and Mbh — cr planes, in a location consistent with 
what observed for pseudobulges. Given the limited growth due to gas accretion, we argue 
that the mass of the central black hole should be above 10^ Mq already at 2 ; ^ 8. The 
effect of AGN feedback on the host galaxy is limited to the very central few hundreds of 
parsecs. Despite being weak, AGN feedback seems to be responsible for the limited growth 
of the central bulge with respect to the original Eris, which results in a significantly flatter 
rotation curve in the inner few kiloparsecs. Moreover, the disk of EiisBH is more prone 
to instabilities, as its bulge is smaller and its disk larger then Eris. As a result, the disk 
of EtisBH undergoes a stronger dynamical evolution relative to Eris and around 2 ; = 0.3 
a weak bar grows into a strong bar of a few disk scale lengths in size. The bar triggers a 
burst of star formation in the inner few hundred parsecs, provides a modest amount of 
new fuel to the central black hole, and causes the bulge of EiisBH to have, by 2 ; = 0, a 
box/peanut morphology. 

Key words: Galaxy: bulge - Galaxy: center - galaxies: active - galaxies: formation - 
quasars: supermassive black holes - methods: numerical 


1 INTRODUCTION 

Present in nearly all spheroids (e.g.. lKormen ^l2004l ). and be¬ 
ing the engines powering active galactic nuclei (AGN), super¬ 
massive black holes (simply black holes, from here onwards) 
seem to be an integral component of massive spheroidal galax- 
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ies. The tight relations between black hole mass and several 
properties of the host spheroid (e.g.. [Ma gorrian et al. 1 9^: 


n ost spn_ . ^ ^ _ _ 

Merritt &: Ferrares3l200l] : lTremaine et al.ll2002l : lHaring &: Rixl 

2 OO 4 I ) give further hints that the life of black holes and their 
hosts are closely linked. Startin g from the first analytical work s 
on black hole “self-regulation” ( Silk &: ReeslIlQQSl : King|l200^ ) . 
a considerable amount of effort has been spent in the last two 
decades in understanding the role of black hole feedback in 
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shaping the properties of the host galaxy and controlling its 
own growth. Significant advance has been made by includ¬ 
ing the physics of black hole accretion and feedback (even if 
limited to sub-grid models) into simulatio ns o f galaxy evolu¬ 
tion. Fo r example, IPi Matteo et al.l (|2005l ) and ISoringel et al.l 
(|2005al ) studied the effects of black hole (thermal) feedback 
during the merger of blue spiral galaxies, and found that 
AGN feedback has a primary role in the color transformation 
of the host galaxy, by quenching star formation and causing 
the elliptical merger remnant to redden in a short timescale. 
The combination of these results and the evidence of scal¬ 
ing relations between black hole mass and bulge properties, 
hint at the idea that is during violent, major merger events 
when black holes acquire most of thei r mass while elliptical 
galaxies form fe.g.. lHoDkins et aI1l2006l ). On larger scales, me¬ 
chanical AGN feedback in the form of jets an d bubbles seems 
to halt cooling flows in galaxy clusters fe.g.. iGhurazov et al.l 
I 2 OO 2 I : ISiiacki SDring^l2006l ) . While black hole growth and 
the effects of AGN feedback have been extensively explored 
in the case of mas s ive early-type galaxies in clusters (e.g ., 
iMartizzi et al.|[2012l : IPubois et al.ll201^ : iMartizzi et al.|[2014 ), 
the role black holes have played in the evolution of blue late- 
type galaxies evolving in smaller haloes is still largely unex¬ 
plored. The small interest the community has devoted to the 
interplay between black holes and spiral galaxies is primarily 
due to the fact that late-type galaxies have small bulges and 
small black holes. Also, their active star formation suggests 
that no quenching is present in those galaxies, and thus AGN 
feedback is not significant. 


In the effort of simulating the evolution of disk galaxies, a 
lot of attention has been instead given to stellar feedback. Ob¬ 
taining realistic late-type galaxies has been a historical chal¬ 
lenge for computational as trophysics since the first simulations 
run in the early nineties. iNavarro BenzI (|l99ll ) found that 
the baryonic component of galaxies evolving in hierarchically 
growing haloes was not able to retain enough angular mo¬ 
mentum to feature a cold disk and flat rotation curves as in 
realistic spirals. Those authors suggested that a proper treat¬ 
ment of supernova feedback (or some other form of heating) at 
early times would be necessary to prevent gas to lose angular 
momentum and catastrophically cool during galaxy interac¬ 
tions. Since these first simulations, our understanding of the 
formation of spiral galaxy has evolved significantly and, in the 
last few years, a large number of works from different groups 
has shown that proper resolution and an accurate treatment 
of star formation and stellar feedback are th e key ingredi¬ 
ents for creating realistic late-type galaxies (e . g., Guedes_^l_al] 
2011 : iBr ook et aT ] I 2 OIII: Aumer et al. l201^: Tokamat3 2013|: 

3 iMari] 


Roskar e t al.ll2o 3 iMarinacci et ^ 2014 : lAgertz Kravtso^ 

2 OI 5 I : iMurante et al.ll2015l ). For further details on thi s topic . 


we refer the reader to the early re view of Maver et al.l , 

the large code-comparison work of lScannaDiecol (|2012l ) and the 
su mmary of all most re cent developments given in section 2 
of Murante et al.l (l2015l ). We discuss here only the work of 
iGuedes et al.l (|201ll ). as that is the starting point of the new 
simulation presented in this paper. iGuedes et al ] (1201 ih per¬ 
formed a zoom-in cosmological hydro simulation of a Milky 
Way-size halo with a rather quiet merger history. They found 
that the combination of high spatial resolution (gravitational 
softening of 120 pc) and the use of a high-density thres hold for 
star formation (as first use d bv |Governa to et al.ll2010l ) and of 
blast-wave feedback (as in IStinson et al.l 20061 ) . contribute to 
the development of a clumpy and inhomogeneous interstellar 


medium, where overlapping SN explosions are able to inject 
enough energy to remove low-angula r momentum material. 
The simulation of IGuedes et al.l (|201ll ). dubbed “Eris”, is one 
of the first successful efforts to produce a realistic late-type 
spiral in a cosmological simulation. 

With the new simulation that we present in this paper, 
“ErisRif ”, we want to examine directly the role of AGN feed¬ 
back in the evolution of late-type galaxies and explore what 
could be the origin and cosmological evolution of the black 
holes that these galaxies host. We do so by adding sub-grid 
physics for the growth and feedback of supermassive black 
holes to the physics already simulated in Eris. The paper is 
organized as follows. Section offers a brief summary of the 
Eris simulation and describes the technical aspects of EiisBH 
. In Section [3] we show the evolution of the black holes in the 
simulation, focusing in particular on the central black hole 
of the main galaxy of the simulation. In Section 0] we show 
how, and to which extent, AGN feedback influences the host 
galaxy, by making a direct comparison between the physical 
properties of EiisBH and those of Eris. Einally, in Section [5] 
we summarize of our results. 


2 THE SIMULATION 

We describe here the technical aspects of the simulation. The 
initial conditions and the physical processes included in the 
simulation are the same as the ones of Eris, which are de¬ 
scribed in Section nm In addition, this new simulation also 
includes assumptions for the formation of massive black holes 
and follows their growth as described in Section [2]2] 


2.1 Eris 

Eris is a cosmological zoom-in A-body/smooth particle hy¬ 
drodynamic (SPH) simulation that follows the evolution of 
a Milky Way-size halo from z = 90 down to z = 0, in a 
Wilkinson Microwav e Anisotropy Probe three-year cosmology 
([Spergel et al.l l2007l ) , with a flat universe with Qm = 0.24, 
Qh = 0.042, ho = 73 km Mpc“^, n = 1 and ag = 0.76. 
It was run with the parallel, sp atially and teniporall y adap- 
tive, treeSPH-code GASOLINE (IWadslev et aljl2004ll for 1.5 
million cpu hours. The halo was selected from the z — 0 out¬ 
put of a low-resolution and dark matter-only simulation of a 
(90 Mpc)^ volume; it was chosen to have approximately the es¬ 
timated size of the Milky Way halcE] and to have had a rather 
quiet merger history, with no major mergers (above the ra¬ 
tio of 1 : 10) after z = 3. The dark matter and the (initial) 
gas particle masses of the high-resolution region of Eris are 
muM = 9.8 X 10^ M© and rugas = 2 x 10^ M©. The gravita¬ 
tional softening length was fixed to 120 physical pc from z = 9 
to the present, and evolved as 1/(1+z) from z = 9 to the start- 
ing redshift (z = 90 ). In brief (for further technical details, see 
IGuedes et al.l 1201 ll ) , the simulation includes Gompton cool¬ 
ing, atomic cooling, metallicity-dependent radiative cooling at 
low temperatures and uniform UV background. Star forma¬ 
tion and stellar feedback are regulated by the star formation 
threshold nsF (set to a value of 5 atom/cc), the star formation 


^ The virial mass if Eris at 2 ; = 0 is Mvir = 7.8 x 10^^ M©, and the 
virial radius R^ \ r = 234 kpc, as given by t he halo finder AMIGA 
(|Gill et al.ll2004l : iKnollmann Sz KnebellSoO^i . 
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Figure 1. Gas density maps in the location where new black hole seeds are about to be inserted in the simulation. 


efficiency esF (set to esF = 0.1) and the fraction of of super¬ 
nova energy that couples to the interstellar medium esN (set 
to 0.8). When the conditions of the gas become favourable 
to star formation (according to t he Schmidt la w), new star 
particles are created following the iKrou 3 (|200lh initial mass 
function, with a starting mass = 6 x 10^ M© (the gas par¬ 
ticle from which the star is created has its mass reduced by 
the same amount). Supernova explosions imply a deposition 
of metals and energy (esN x 10^^ erg) to the nearest neigh¬ 
bour gas particles. The affected gas has its co oling shut off 
until t he end of the supernova blastwave, as in IStinson et al.l 
(I 2 OO 6 I). Wind s and stellar mass loss are also modelled as in 
IStinson et al.l (|2006h . 

As discussed in (|Guedes et al.|[201ll T a high value for the 
density threshold for star formation parameter was possible 
thanks to the high mass and spatial resolutions of the simu¬ 
lation, which allow to resolve the clouds where star formation 
occurs. With this choice, star formation takes place in confined 
regions, giving rise to a clumpy, inhomogeneous interstellar 
medium, where overlapping supernova explosions inject en¬ 
ergy in a localized manner. This localized energy injection is 
able to create galactic outflows which expel low-angular mo¬ 
mentum material. 

The formation of a clumpy and inhomogeneous medium 
seems to be the key in gredient to form a r ealistic late-type 
spiral. As described in (|Guedes et al.l [ 2 OI 1 I 1 . in fact, Eris is 
consistent with a large number of observational aspects of the 
Milky Way, from its structural properties to the mass content 
of its different components. Just to give few examples, it has 
a fiat rotation curve, a low photometric bulge-to-disc (B/D) 
ratio, it falls on the Tully-Fisher and stellar-mass/halo-mass 
relation and it has a baryonic mass fraction within the virial 
radius which is 30% lower than the cosmic value. 


2.2 EivisBH : Eris plus black hole formation, growth 
and feedback 

ErisBH is a replica of Eris, but with additional prescriptions 
for the seeding and growth of supermassive black holes, and 
thermal feedback during gas accretion. 


2.2.1 Seeding proeedure 

The origin of black hole seeds is still largely a mystery and sub- 
je ct of intense theoretical inv estigation (see, e.g., the review 
of IVolonteri Bellovarvir2012l l. Thus, the criteria for inserting 
black h ole seeds in a cosmologic al si mulation is somewhat ar - 
bitrary. IPi Matteo et al.l (|200^ 1 and iBooth Schavd (|2009l l , 
for example, have chosen to insert a seed black hole in every 
dark matter ha l o that rises above a given threshold in mass. 
iBellovarv et al.l (|2010l T instead, chose to connect the forma¬ 
tion of black holes to star formation: gas particles with zero 
metallicity and with density above the set threshold for star 
formation, have a certain probability to be transformed into 
black holes rather than stars, where this probability parameter 
is tuned to reproduce the black hole seed halo occupat ion prob¬ 
ability at z = 3 suggested by IVolonteri et al.l (|2008l l (which, 
however, is based on a very specific formation model for black 
hole seeds). Here, we impose both a requirement on resolution 
and high-density gas environment for the introduction of new 
black holes. A seed is inserted in all systems that satisfy the 
following conditions: 

• are bound, according to the AMIG A halo finder 
(|Gill et al.l l200Jj : iKnollmann Knebd l2009l ). and resolved 
with at least 10^ particles, 

• have at minimum of 10 gas particles with density above 
100 atoms/cc. 
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z - seed 

Ngas > 100 atom/cc 

Mvir 

Mgas 

Afstar 

MBHseed 

1 

8.5 

1354 

7.2 X 10^ M© 

7.2 X 10® M© 

1.7 X 10® M© 

8.7 X 10^ M© 

2 

5.8 

579 

4.8 X 109 Mq 

5.9 X 10® M© 

1.7 X 10® M© 

1.3 X 10® M© 

3 

4.7 

318 

7.0 X 10® Mq 

1.4 X 10^ M© 

2.0 X 10® M© 

7.6 X 10® M© 

4 

3.9 

13 

6.4 X 10® Mq 

7.0 X 10® M© 

1.2 X 10® M© 

0.8 X 10® M© 


Table 1. Summary of the physical properties of the systems that qualify as proper sites for black hole seeding, following the conditions 
described in Section mm The first column gives the redshift at which the black hole is inserted in the system. The second column gives 
the number of gas particles that have a density higher than 100 atom/cc. The remaining columns give the virial, gas and stellar mass of the 
systems as well as the mass of the black hole seed. 


The resolution criterion is motivated by past works on numeri¬ 
cal convergence of angular momentum transport in astrophys- 
ical disks simulated by SPH, which have shown that at least 
10^ particles are needed to keep a number of numerical effects 
under control, such as spurious hydrodynamical drags as well 
as enhanced gravitation al torques from a noisy halo potential 
(|Kaufmann et al]l20Q7l l. 

Galaxies that satisfy the conditions listed above are then 
seeded with a black hole, provided that they do not host one 
already. This is done by selecting the star particle that is clos¬ 
est to the high-density gas particle with the deepest potential 
and converting it into a sink particle. The initial black hole 
mass is set to be proportional to the number of high-density 
gas particle^, so that the higher is the density of the gas, and 
the larger is the high-density region, the more massive is the 
newly-formed black hole. 

From the beginning of the simulation down to ~ 3, only 
six structures raise above the mass (or number of particles) 
limit that we have imposed. Of those, four satisfy also the 
condition on the gas density. At lower redshift, the gas densi¬ 
ties become generally lower, and we stop looking for possible 
sites for new black holes seeds. 

The first system which satisfies the seeding conditions is, 
at z ~ 8.5, the progenitor of the main galaxy of the simula¬ 
tion. The other three seeds are inserted in satellite galaxies at 
lower redshifts. Figure [T] shows maps of the gas density in the 
regions where new seeds are about to be added to the simu¬ 
lation. Table [T] gives further information on the properties of 
the systems where black holes are inserted. Given the imposed 
limit on the number of particles that a protogalaxy needs to 
have to host a black hole, the masses in gas, star and dark 
matter do not differ substantially in the four galaxies at the 
time each black hole is inserted. Seed masses (last column in 
the table) also do not differ significantly, with about one order 
of magnitude separating the most to the least massive seed. 
As a consequence, all black hole seeds are located close to each 
other on the Mbh— stellar mass plane, slightly above the ex¬ 
trapolation to lower masses of the relation established by the 
most massive black holes (see left panel of Figure [T]). 


2.2.2 Black hole accretion and feedback model 


After being inserted into the simulation by converting the ap¬ 
propriate star particle into a sink-type particle, black holes 
start accreting gas isotropically from the surrounding medium, 
following the widely used Bondi-Hoyle-Lyttleton formalis m 
(|Hovle Lvttlet^ll939l : iBondi Hovlel[l944l : lBondil[T952h : 


Afsondi — 


AirG^Mi^p 

(ci+u2)3/2’ 


( 1 ) 


where p and Cg are, respectively, the density and sound speed 
of the gas, Mbh is the mass of the black hole, and v is the 
velocity of the black hole relative to the gas. Growth pro¬ 
ceeds according to Equation [U but the maximum allowed ac¬ 
cretion rate is set to the Eddington accretion rate MEdd = 
{^TrGMBHmp) f{r]aTc) (where nip is the proton mass, c is the 
speed of light, ctt the Thomson cross-section and rj is the ac¬ 
cretion efficiency, assumed to be 0.1), so that: 


Mbh 


-^i^Bondi if -^i^Bondi ^ -^i^Edd 
MEdd if MBondi > MEdd 


( 2 ) 


Eeedback from the accreting black hole is modeled by assum¬ 
ing that a fraction e/ = 0.05 of the radiated luminosit>0 is 
converted into thermal energy th at heats the gas surroun ding 
the black hole . Initial ly used bv IPi Matteo et "all (|2005l i and 
ISoringel et al.l (|2005bl i in their study of the effects of black hole 
growth and feedback in isolated galaxy mergers, the Bondi- 
prescription for estimating the growth of massive black hole 
has been used widely in the community, also in large cosmo¬ 
logica l volumes (e.g. IPi Matteo et alJ l2008l : iBooth fc Schavj 
l2009l^ (we discuss the limitation of the Bondi-Hoyle-Lyt tleton 
formalism below in section 13.2D . iBellovarv et al.l (|2010h have 
included this prescription into the GASOLINE code to study 
wandering black holes, that is, black holes that are the rem - 
nants of stripped satellite cores. As in IBellovarv et all (|2010I L 
black holes are allowed to merge if they are within one anothers 
softening length and if they fulfill the criterion | < AaAr, 

where Av and Aa are the differences in velocity and acceler¬ 
ation of the two black holes, and Ar is the distance between 
them. 


^ specifically, we calculate the number Ngas,dens of gets particles 
that have a density larger than 50% the density of the most dense 
gas particle in the galaxy. The mass of the black hole seed is then 
given by Ngas,dens times the mass of the converted star particle. 


3 The bolometric radiated luminosity is given by LboI = 
^radMBHC^, where Crad is the radiative efficiency. For simplicity, we 
assume here e^ad = Vi which is a good approximation at high accre¬ 
tion rates when black holes are likely g rowing from geometrically - 
thin and optically-thick accretion disks (IShakura &: SunvaevlllOT^f . 
but it might be an overestimate of the radiative power when 
accretion rates are hi ghly sub-Eddington (IChurazov et al.l l2005l : 
iMerloni Sz Hein'3l2008lV 
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Figure 2. Upper panel: mass evolution of the four black holes in the 
simulation. Symbols stop when the corresponding black hole merges 
with the black hole of the central galaxy (yellow track). Lower panel: 
again mass evolution of the four black holes, but now normalized by 
the respective initial mass. The vertical arrows indicate the redshift 
at which black hole mergers take place. 


3 BLACK HOLE PROPERTIES 
3.1 Black hole growth 

The upper panel of Figure [2] shows the mass evolution of the 
four black holes in the simulation. Curves stop when the cor¬ 
responding black hole merges with the black hole of the cen¬ 
tral galaxy (yellow symbols). The mergers (marked by the ar¬ 
rows) are clearly visible as “jumps” in the yellow track. The 
lower panel of the figure shows again the growth of the four 
black holes, this time normalized by the corresponding ini¬ 
tial mass. Growth by gas accretion is clearly limited in all 
black holes. The central black approximately doubles its mass 
though mergers and, since the time it was seeded into the sim¬ 
ulation at z ~ 8.5, it grows, in total, only by a factor of ~ 3, 
reaching a final mass of about 2.6 million Solar masses. The gas 
growth we find in our simula tion is much more mod est com¬ 
pared to the one obtained bv iMarinacci et ^ (|2014l ) in their 
simulations of eight Milky Way-sized haloes: all black holes 
hosted by their simulated central galaxies reach final masses 
around 10® M©. We argue that the differences in the results 
are likely primarily due t o the different resolution s adopted in 
our and their simulations: iMarinacci et~^ li^ have a grav¬ 
itational softening of about 700 pc, while we resolve down to 
120 pc. We argue that the lower resolution of their simulations 
implies an overestimate of the gas supply available for accre¬ 
tion. In our simulation, for example, the amount of gas inside 
~ 1 kpc is around 10® M©, almost an order of magnitude 
higher than in the inner 200 — 300 pc. 

Figure [3] shows the gas accretion rate evolution of the cen¬ 
tral black hole. The orange points in the upper panel indicate 
the instantaneous accretion rate, while the orange bullets show 
the median over several timesteps (the error bars indicate the 
16 and 84 percentiles of the distribution). The median accre¬ 
tion rates decrease only gently with redshift, and are typically 
quite low, ~ 10“^ — 10“^ M©/ yr, but there are periodic fluc¬ 
tuations and peaks at higher accretion rates. While there is a 
global gentle decrease in the median accretion rates with de¬ 
creasing redshift, at 2 : ~ 0.3 the median rate increases again 
of about an order of magnitude: as we will discuss in the next 
section, this is approximately the redshift at which the disk 
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Figure 3. Upper panel: evolution of the instantaneous accretion 
rate of the black hole in the main galaxy (orange points) with me¬ 
dian values over multiple timesteps (orange bullets) and 16 and 84 
percentiles (error bars). The red curve indicates what would be the 
Eddington rate for the black hole at each time. The orange dia¬ 
monds highlight when the Bondi formula would have given a super- 
Eddington rate, given the gas properties around the black hole. The 
grey squares show the evolution of the sfr (multiplied by 10“^) in 
the inner 2 kpc of the galaxy. Lower panel: evolution of the Edding¬ 
ton ratio, with again median values over many timesteps and 16 
and 84 percentiles (filled bullets with error bars). The dashed line 
indicates the Eddington limit and, again, the diamonds show when 
the accretion rates would have been super-Eddington. 



of the galaxy is experiencing an instability event that leads to 
the formation of a strong bar visible in the galaxy at z = 0. 
During the process of bar formation, the gas mass and the star 
formation increase in the very central few hundred parsecs and 
the conditions around the black hole also become favourable 
for an increase of the accretion rate. By z = 0, the accretion 
rate has lowered again to values around 10“^ M©/ yr. 

In the same Figure, the red bullets indicate, for reference, 
which would be the instantaneous Eddington accretion rate for 
the black hole. We also explicitly show when the conditions of 
the gas around the black hole would have given a Bondi accre¬ 
tion rate higher than the Eddington limit (orange diamonds): 
this happens rarely, and only at high redshift. However, as we 
do not allow super-Eddington accretion, the accretion rate in 
those rare cases is set to the Eddington limit. 

The bottom panel explicitly shows the evolution of the 
Eddington ratio {fEdd = MbhIME dd)^ which is typically be¬ 
tween 10“^ and 10“^. Again, the diamonds show the few rare 
times where the Bondi formula would give accretion rates 
higher than the Eddington limit. By z = 0 the Eddington 
ratio reaches the lowest median values of ~ 10“^. 

Accretion rates between 10“^ and 10“® M©/ yr are, con¬ 
servatively, the accretion rates expected from Bondi accretion 
of hot gas and stell ar mass los s in the center of nearby galax¬ 
ies, as discussed by 113 (l200flh . At z = 0 our simulated black 
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hole is accreting close to those values, at ~ 10“^ M©/ yr, 
which, assuming a radiative efficiency Crad = 0.1, corresponds 
to a bolometric luminosity of ~ 10^ L©. This is close to the 
median v alue of nu clear luminosity of local Seyfert galaxies 
found by I Hoi (|2009l b As discussed in Section [2.2.21 however, 
the radiative efficiency is likely lower than the 10% value we 
have adopted here, and can be as low as Crad ~ 10“^ — 10~^ 
for the accretion rates we have in our simulation at late time. 
At those low rates, in fact, energy from the accreting black 
hole is likely released a s kinetic, rather than r adiative, power 
(jChurazov et al.l [20051 : iMerloni Heinjl2008l b With signifi¬ 
cantly lower efficiencies of energy release, the luminosity of the 
nucleus of ErisBH could be several order of magnitu de lower, 
thus being closer to the luminosities that IHoI (( 2 OO 9 I I find for 
a large fraction of local galaxies (including the Milky Wa^Q). 
Had we assumed a lower value for the radiative efficiency for 
the phases at the lowest Eddinton ratios, the effect of radia¬ 
tive AGN feedback on the galaxy would have been even lower 
than what we have obtained in our simulation. 

In the top panel of Figure [3] we also show the evolution of 
the star formation rate (multiplied by 10“^ for plotting pur¬ 
poses, and with values averaged over 10 snapshots, which cor¬ 
respond to about 300 Myr) in the inner 2 kpc of the galaxy. 
As the bulge stars dominate approximately the inner 2 kpc 
of the galaxy at = 0 (see section JilL a nd as most bulge 
stars are formed in-situ fsee iGuedes et al1l2013l b this sfr is 
a good proxy for the sfr of the stars within the bulge. The 
black hole accretion rate is between 6 and 4 orders of magni¬ 
tude smaller than the star formation rate. As expected for a 
late-type galaxy, the growth rate of the central black hole is 
thus significantly smaller than the growth rate of the bulge. 
In section [ 33 ] we explicitly show how this translate into black 
hole-galaxy scaling relations. 

Given the little gas supply, the value of the initial seed 
mass is rather important. Generalizing our results, we expect 
the progenitors of Milky Way-size galaxies to be hosting an 

8. For seeds from 


PopIII remnants (e.g..lBond et al.lll984j 

iMadau & Reesll200ll: 

iHaiman & Huill200ll: iTanaka & Haimanl 

l2009lb edven the lit- 


tie time available before z 8, uninterrupted growth close 
to the Eddington rate would be required to reach interme¬ 
diate masses (unl ess super-Eddington accretion is possible , 
see, for example, Wvithe Loe^ 20121 : iMadau et al.l l2014l : 
IVolonteri et al.ll2015l b An alternative are seeds from direct- 
collapse, which could be already between 10^ and 10® M© 
at the time of formation. Given the properties of Milky Way 
progenitors, dire ct collapse black holes born in metal-free pro¬ 
togalaxies (e.g., iLodato Nataraian I2OO6I: IWise et al.l I2OO8I: 


Regan Haehne ltl l2009l: Johnson et al.l 2011 : Agarwal et al l 
2 OI 2 I : IPiikstra et al.ll2014l l are a more plausible scenario than 


the direct collapse black holes formed after galaxy major merg¬ 
ers, as these events are only possible during mergers of much 
more massive gala xies (jMaver et al.lI 2 OIOI : iBonoli et al.ll2014l : 
iMaver et al.|[20l3 b 


The estimated lu minosity of the Gala c tic center is L^g] < 
10^^ erg/s (see,e.g . I Skinner et al.l Il987l : [P^linsk^^^^ ll99j; 
iBaganoff et aklBoO^l and other references in lNamvar^^iLrjlQQSl V 


3.2 A discussion on the Bondi-Hoyle-Lyttleton 
prescription 


As introduced in the methodology section, in this paper we 
have used a specific model for black hole accretion, the Bondi- 
Hoyle-Littleton model, which is simple and widely used in 
the literature of cosmological simulations. At the same time, 
though, it is known to be often inaccurate for flows in the nu¬ 
clei of galactic disks, which are neither spherical and hydrody- 
namical, but are rather governed by the effects of gravitational 
torques in the re distribution of angular momentum and ki¬ 
netic energy fe.g. JPebattista et al.ll2006l:lHoD kins fc QuataertI 
I 2 OIOI . l201lh . In particular. iHopkins QuataertI (|201lh have 
proposed a different sub-grid model for accretion onto mas¬ 
sive black holes in simulations that is based on the analysis 
of angular momentum transport in galactic disks of stars and 
gas, where the self-gravitating perturbation destabilizing the 
axisymmetric gas flow is assumed to be coming from the stellar 
potential only. The latter is a good approximation in our case 
since star formation is vigorous in the nucleus at high redshift 
(see Figure [T^ . leading to stellar-to-gaseous mass ratio larger 
than 10 at distances below 300 pc from the center already 
at redshift z 4 (see Figure [8l). iHopkins QuataertI (| 20 lib 
have shown that their proposed new sub-grid model captures 
very closely the mass inflow rates at small scales occurring in 
very high resolution simulations, which resolve down to ~ pc 
scales, hence close to the boundary of the accretion disc, while 
the Bondi model can overestimate by even more than an order 
of magnitude the inflow rate occurring in simulations at such 
scales. Since we only resolve gravity to about 100 pc scales, 
the size of our gravitational softening length, we can apply 
their sub-grid model at such scales. According to this model, 
the mass inflow rate M onto the black hole is ~ aMgas/Torb, 
where Mgas is the gas mass at the smallest scale, Torb is the 
orbital time at the same distance and a is a constant related 
to the amplitude of gravitational torques, which is maximum 
in major mergers (a = 1) and can be as small as a = 0.01 in 
weakly self-gravitating disk. We can apply the model to the 
high-redshift phase of black hole growth (z = 2 — 4) which is 
the one yielding, on average, the highest accretion rates (see 
Figure [3|). Assuming a a = 0.01 — 0.1, since most of the time 
our galaxy does not undergo significant mergers, and taking 
Mgas a few times 10^ M© (corresponding to the values in the 
original Eris simulation shown in Figure (S]) and Torb ~ 10^ yr, 
a typical value for the inner few hundred parsecs, we obtain 
M ~ 10~^ —10~^M© yr. Those values for the accretion rate are 
higher than the accretion rates measured in the ErisTif simu¬ 
lation, despite it adopts the Bondi model. This is reassuring, 
since feedback is expected to decrease the accretion rates fur¬ 
ther, yet there is no evidence that the Bondi model employed 
is boosting severely the accretion rates. This suggests that in 
our case the accretion rates are not strongly dependent on the 
accretion model adopted, but they are rather bound to be low 
because they arise from the weak gas inflows occurring at the 
resolved scales, from few kpc to a few softening lengths. This 
is line with the notion that we are modeling a fairly quiescent 
halo that does not undergo many prominent mergers, and this 
seems to naturally lead to the formation of a late-type spiral 
with a pseudobulge. It also confirms that our general key con¬ 
clusion that massive black holes in late-type spirals grow little 
by accretion should be quite general. 


© 0000 RAS, MNRAS 000 , 000-000 


































































Black Hole Evolution in a Late-type Galaxy 7 



Figure 4. Time sequence of the last black hole merger in the simulation. Black holes are represented by the yellow and red bullets. The gas 
map is color-coded according to the particle density. In black are overplotted stellar isodensity contours. 
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Figure 5. Details of the spatial evolution of the two black holes 
shown in the maps of Figure 0] Upper left panel: evolution of the 
x—y position of the two black holes. Top right panel: x—y projection 
of the orbit of the secondary black hole around the primary. Lower 
panel: redshift evolution of the separation between the two black 
holes. 


3.3 A close look to black hole mergers 

As previously discussed, mergers with black holes from satel¬ 
lite galaxies contribute significantly to the mass-growth bud¬ 
get of the black hole hosted by the central, most massive 
galaxy of the simulation. The resolution of the simulation al¬ 
lows us to follow the evolution of black hole binaries down to 


few hundreds of pc. The maps in Figure |3]show the evolution 
of the last black hole merger in the simulation, occurring at 
z ^ 1, with gas particles colored according to their density 
and black lines indicating stellar isodensity contours. At the 
beginning of the sequence, the stellar core hosting the satel¬ 
lite black hole is still visible in the isodensity contours. When 
the remnant of the satellite galaxy is finally disrupted, the 
secondary black hole is dragged to the center of the primary 
galaxy to eventually merge with the more massive black hole. 
This sequence spans ~ 250 Myr in the cosmology assumed 
here. 

Figure [5] shows quantitatively the position, separation and 
orbit evolution of the two black holes from several kpc scales 
down to the softening scale of the simulation. The upper left 
panel shows the x — y position of the two black holes in the 
same time interval shown in Figure 0] where the positions have 
been centered at the location of the primary black hole at the 
time we start tracking the binary system. While the central 
black hole (orange curve) wiggles around its original position, 
the satellite black hole shows a typical spiraling pattern until it 
finally joins the central black hole. The top-right panel shows 
instead a x — y projection of the orbit of the satellite black 
hole around the primary one and, finally, the bottom panel of 
the same figure shows the redshift evolution of the separation 
between the two black holes. In about 250 Myr the black holes 
go from a distance of several kpc to a separation smaller than 
the spatial resolution of the simulation, at which point they 
are assumed to merge. Such short timescale for the evolution 
of the black hole binary system is consistent with the results of 
high-resolution mer ger simulations of, e.g., (|Maver et al.ll200'^ : 
IChapon et ahllioisl i. 

However, we also note that detailed numerical simulations 
of massive black hole mergers in circumnuclear disks have 
shown that the formation of a tight binary below our reso¬ 
lution scale of 120 pc can be del ayed to up to 10^ y r in the 
favourable case of major mergers ((Roskar et al. I [20151 ). and to 
even longer timescales (up to a Gyr) in the case of minor merg¬ 
ers, due to a variety of processes such as ram pressure strip- 
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Figure 6. Mbh ~ cr plane: th e black symbols a r e obs ervational 
points from the compilation of iMcConnell Sz mU (|2013l) (best fit 
shown by the dotted black line), the orange diamond is the location 
of the Milky Way in the plane while the yellow circle indicates the 
location of KrisBH . 


ping of the secondary galaxy co r e surr ounding the secondary 
black hole (jCallegari et al.|[2009l . l201lh . In any case, all black 
hole mergers in our simulation occur at very early times, when 
dynamical timescales are still very short, possibly explaining 
that neglecting to capture the small scale dynamics is not an 
issue. We finally note that, with our resolution, we encounter 
no problems in following the black hole orbital decay directly 
down to the scale of the softening leng th, namely without h av- 
ing to impose an artificial drag force (|Tremmel et al.l[ 20 1 5l i . 


3.4 Black hole-galaxy scaling relations 

Since their dis cov ery m ore than 20 years a go (e.g. , 
Magorrian et al.l _ 19981 : iMerritt Ferrari I 2 OOII : 


Tremaine et al.l 12002 )^ the tight relations between black 


hole mass and various properties of the host galaxies have 
been interpreted as an evidence of some type of co-evolution 
between massive black holes and their hosts. Those relations 
have been first defined for massive bulges or elliptical galaxies, 
and has not been clear whether they hold when going to low 
mass galaxies or galaxies with di f ferent morphologies. Recent 


studies (e.g., Greene et al.l I 2 OO 8 I: Hu 


2011 


Ijiang et~al]l201ll: iMathur et al 


2012 


200g; [Kormend^j^^ 


iGraham Scott 


20151 ) have found that black holes in pseudobulges do not 


correlate in the same way, if at all, with galaxy properties 
as black holes in class i cal b ulges. Studying active low-mass 
galaxies, iGreene et al.l (|2008l) found, for example, that black 


holes in pseudobulges sit on the extrapolation to lower masses 
of the Mbh — cr relation, while they lie almost an order of 
magnitude below the Mbh — Mnuige relation. The different 
behaviour for classical bulges and pseudobulges seems to be 
due to the different relation be tween stellar mass and velocit y 
dispersion that t hey follow (|Gadotti KauffmannI l2009l ). 
iKormendv et al ] il 2 niih studied systems with dynamically- 
estimated black hole masses, and found evidence that no 
correlation actually exists between black hole masses and 
pseudobulges. Those authors argue that this is likely due 
to the different history of black holes evolving in pseudob¬ 


ulges and the ones in classical bulges: black hole growth in 
pseudobulges is likely driven by small-scales and stochastic 
events and is typically highly sub-Eddington, in contrast 
with the growth of the most massive black holes, which is 
driven by large scale dramatic events, such as the violent 
mergers, that are able to both efficiently feed black holes and 
to generate classical bulges and elliptical galaxies. The bulge 
of EtisBH has the properties of a pseudobulge rather than 
a classical bulge (see Section oil . similarly to the bulge of 
Er is, whose evolution has been studied in a dedicated paper 
by iGuedes et al.l (|2013l ) . We can thus directly see whether our 
black hole has properties consistent with the ones of black 
holes in observed pseudobulges. 


In Eigure[6]is shown where the EiisBH central black hole 
sits in the Mbh — cr relation (yellow bullet ), compared with 
the ob servational data points compiled by iMcGonnell &: Mai 
(|2013l ). We estimated the velocity dispersion for EvisBH us¬ 
ing the I/V 3 of the 3 — D velocity dispersion within 2 kpc 
from the center. With this assumption, we obtain a value for 
cr of 109 km s“^. Our calculation for the velocity dispersion 
gives values that are quite comparable to obs ervational long- 
slit m easurements of face-on nearby galaxies (jBellovarv et al.l 
I 2 OI 4 I ). The obtained value of 109 km s“^ is also within the 
errors on the velocity dispersion of 103( zb20) km s~^ for 
the M ilky Way given in the compilation of lMcGonnell &: Mai 
(|2013l ). We see that ErisHifsits slightly below the relation, as 
its black hole mass is about 2/3 the mass of t he Milky Way cen- 
tral black hole, assuming th at SgrA*, as inlMcGonneh Ma 

S ), is 4.1 X 10® Mq (from: [Ghez et al.ll2008l : lGillessen et al. 

. 


Eigure[3shows instead the scaling relations between black 
hole mass and stellar mass and between black hole mass and 
the bulge component of galaxies. He re we use as reference ob - 
servational points the compilation of lErwin' Gadottil (|2012l ) 
(red and blue symbols, representing the location of early and 
late type galaxies respectively). As in the previous hgure, the 
orange diamond represent the location of the Milky Way and 
the yellow bullet the location of EiisBH . In Section 14.21 we 
perform a bulge and disc decomposition by htting a double- 
Sersic prohle to the surface density prohle of our simulated 
galaxy: integrating the best-ht parameters of the Sersic pro- 
hles, we obtain estimates for the bulge and disk masses (see 
Table 14.21) . The stellar mass used here is the sum of the bulge 
and disc derived from the decomposition. EiisBH sits well be¬ 
low the relation, even though the black hole seeds start off 
close to the extrapolation of the relation at lower masses: this 
is because the black hole accretion rates are much lower than 
what would be required to keep up with the star formation 
rate in the central region (see Eigure[3|). 


Having only one simulation, we can not give quantitative 
predictions on whether black holes in pseudobulges follow a 
different relation or do not follow any relation at ah with the 
properties of the host. But we note that the consistency with 
the Mbh — cr relation on one side, and the deviation from the 
Mbh — Mn uige on the other sid e, are consistent with what 
is found by Greene et H] ll2008h. but i s also consistent with 
the picture drawn bv iKormendv et all (|201ll ). given the slow 
Seyfert-like growth of our simulated black hole. 
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Figure 7. Mbh — -^star and Mbh — A^Bulge relations. The red and blue symbols are from a compilation of lErwin fc Gadottil 1 I 2 OI 2 I I for 
early type and late type galaxies, respectively (with fits for the Mbh — MBu ig-e rel ation shown by the dashed lines). The orange diamond is 
the location of the Milky Way in the plane as reported bv lErwin Sz Gadottil (l2012l V The yellow circle indicates the location of ErisBiif . The 
colored bullets on the left of the plot show the location of the seeds, color-coded as in Figure ^ 



z 
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Figure 8. Redshift evolution of gas mass (left panel) and star mass (right panel) enclosed within different radii from the galaxy center. The 
filled circles are for EjTisBH , while the open ones for Eris. The bottom panels show the ratios of the corresponding enclosed mass between 
Eris and KrisBH . 


4 BLACK HOLE FEEDBACK AND ITS 
EFFECTS ON THE HOST GALAXY 

As described in Section 12 .2.21 the simulation includes also 
a prescription for black hole feedback: a fraction of the en¬ 
ergy extracted in the accretion process couples thermally with 
the surrounding medium. This standard prescription has been 
used extensively in other works, and its effects led to the idea 
of a “self-regulation” of black hole growth to explain the ori¬ 
gin of the Mbh — cr relation and the ejection of gas from the 
host galaxy and the resulting quenching of star formation (e.g., 


iDi Matteo et al.ll2005ll . Given that the accretion rates are gen¬ 
erally quite small for the central black hole of our simulated 
galaxy, we expect the effects of feedback to be modest. We 
find, in fact, that the amount of energy released by the black 
hole is significantly lower than the one released by supernova 
explosions (integrating from z ~ 3 to 2 : = 0, the total energy 
input from supernovae, both of type la and II, is more than 
an order of magnitude higher than the energy input from the 
central black hole). Despite not playing any dramatic role in 
shaping the host galaxy, AGN feedback seems, however, to be 
able to keep hot the very central region of the galaxy, pre- 
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Figure 9. Gas density maps at z = 0.5 (upper panels) and at 2 : 


0 (lower panels) for ErisBH (left) and Eris (right). 


venting the bulge to grow in the same way as in the original 
Eris. 

To isolate the effects of AGN feedback on the EiisBH 
galaxy, we compare the properties of its gaseous and stellar 
components directly with the Eris galaxy, where AGN feed¬ 
back was ignored in the calculation. 

4.1 Gas component 

In the left panel of Eigure [8] we show the redshift evolution 
of the total mass in gas of Eris and ErisBH within spheres of 
different radii from the galaxy center. At high redshift (above 
z ~ 0.5), the gas content of the two simulations is very similar, 
except in the very inner region (purple curves), where the gas 
budget of EiisBH is systematically lower than the one of Eris, 
a consequence of the weak feedback from the slowly growing 
black hole. 

While at large radii the total mass in gas remains very 
similar in the two simulations down to the present time, in 
the inner 2 kpc of ErisBH (cyan curves) the amount of gas 
starts decreasing at low redshift and, by z = 0, it is about a 
third of the one of Eris. The steep decline in gas content at 
those scales starts around z ~ 0.3 — 0.4, which is when the 
strong bar visible in the galaxy at z = 0 starts forming (see 


next section); as the bar gets stronger, it “clears” the center of 
the galaxy of its gas content. An exception is the very central 
region (the inner few hundred parsecs, purple curves), where 
the amount of gas slightly increases at those late times; this 
raise of gas mass in the center of the galaxy results in a slight 
increase of the black hole accretion rate (see Eigure [S]) as well 
as of the star formation rate (see Eigure [12]). 


In Eigure [9] we show density maps at z = 0.5 (upper 
panels) and z = 0 (lower panels) in a region of 18 kpc of ra¬ 
dius around the center of EiisBH (left column) and Eris (right 
column). The angle of projection is chosen so that the disks 
appear face-on. While, as just discussed, the total amount of 
gas at large scales is the approximately the same in both sim¬ 
ulations, gas seems to be more diffuse in the Erisgalaxy, as 
also its stellar disk is larger than the one of Eris (see the next 
subsection). But the most striking difference between the two 
simulations is probably the “empty” region, of about 2.5 kpc 
of radius, in the center of EiisBH at z = 0. This region has 
approximately the physical extent of the bar that has formed 
at z ~ 0.3 and that, as we mentioned above, by z = 0 has 
cleared the center of the galaxy. 
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Figure 10. Star density maps at 2 : = 0.5 (upper panels) and at 2 : = 0 (lower panels) for ErisBH (left) and Eris (right). 
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Figure 11. Zoom into the central region of ErisBH (left) and Eris (right) at 2 ; = 0, now viewed with the disk edge-on. 
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Figure 12. Comparison between the star formation history of 
ErisBH (filled squares) and Eris (empty squares). The black sym¬ 
bols show the evolution of the total sfr ( Mq/ yr) in the simulation, 
while the colored symbols show the sfr within a given radius from 
the center of the galaxy, as stated in the legend. 
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Figure 13. Age distribution of stars located, at 2 ; = 0, in shells with 
radial distance from the galaxy center as defined in the caption. 
Solid lines are for ErisBH and the dotted ones for Eris. 


4.2 Stellar component 

Figure [10] shows maps of the stellar density projected face- 
on, for the same regions of Figure (9] At z = 0.5 the main 
differences between the two simulations are the slightly larger 
disk and the smaller bulge of ErisBH . It is again at z = 
0 where we see the most remarkable differences between the 
two simulations: in ErisBH the bulge is not only smaller than 
the one of Eris, but now it also features a strong bar, clearly 
visible in the map, and which has a radius of about 3 kpc. The 
presence of the bar causes the bulge to have a boxy-peanut 
morphology when viewed edg e-on, as shown in Fig ure [TTI 

As discussed in details in iGuedes et al.l (|2013 I f the pseu¬ 
dobulge of Eris is the result of an inside-out growth, mainly 
evolved from a stellar bar that formed at high redshift and 
that has progressively weakened after z = 1. In ErisBH the 
evolution at recent times is different: the weak AGN feedback 
has inhibited the growth of the bulge, while its disk has actu¬ 
ally grown to a larger size with respect to Eris (see the right 
panel of [8] for quantitative differences in the stellar content 
of the two simulations at various radii). Being the disk larger 
and the bulge smaller (and thus with less stabilizing power), 
the disk of ErisBH is more prone to instabilities. While the 
bar of Eris thus progressively weakens after z = 1, in ErisBH 
the disk becomes again unstable at z ~ 0. and the bar gains 
again strength with several consequences for the stellar and 
gas properties (e.g., EiguresHKHnSj). 

The stellar density profiles of ErisBH and Eris are shown 
in Eigure (M] together with the best-fits for the bulge and 
disk components. The profiles have been fitted with a double 
Sersic profile, with the Sersic index for the bulge left uncon¬ 
strained, and the one for the disk fixed at 1 (exponential disk). 

^ the ratio r^peak/ \/ / R(i (where r’peak is the peak of the cir¬ 
cular velocity and Mf^andR^ are the mass and scale-length of the 
disk) is tradition ally used to estimate the stability of galactic disks 
i|Mo et ahlflQQSli . Simulations show that a disk becomes unstable 
when the ratio becomes smaller than 1.1. We find that, at 2 : = 0.3, 
ErisBH has precisely a value of 1.1, while, for comparison, Eris 
has a value of 1.4. Those numbers are consistent with the disk of 
ErisBH being about to become unstable, while Eris does not suffer 
any instability event. 


The best values for the fits are given in Table The bulge of 
ErisBH features a smaller Sersic index than Eris (ng = 0.9 and 
Us = 1.2, respectively), with a value lower than unity, typical 
of boxy bulges. The scale radius of the disk of ErisBH is in¬ 
stead larger than the one of Eris (i^disk = 2.9 and i^disk = 2.3, 
respectively). We obtained the total bulge and disk masses by 
integrating the Sersic profiles, and got the values also given 
in Table 14.21 the bulge of ErisBH is almost half the one of 
Eris, while the disk is slightly larger. As discussed above, the 
weak AGN feedback is responsible for the lower reservoir of 
gas in the central few hundred parsecs of ErisBH compared 
to Eris and the consequent lower star formation rate and stel¬ 
lar content in the central region of the galaxy (see Eigures 
[T^ and 0). We note that the values we obtained here for the 
Sersic index of the bulge and the scale radii for Eris differ from 
the ones given in (jGuedes et al.|[201ll i. There, the fit was per¬ 
formed with the GALEIT code run on the surface brightness 
profile obtained after processing th e simulated galaxy with 
SUNRISE. Using this methodology, iGuedes et al.l (|201lh ob¬ 
tained a photometric B/D ratio of 0.35 in the B-band, lower 
of about a factor of 2 than the B/D ratio we obtain here. 
Assuming that for ErisBH the ratio between the photometric 
B/D ratio and the B/D ratio as we calculated here is the 
same as for Eris (which is plausible, given the similar stellar 
ages at galaxy scales between the two simulations, see Eigure 
dsi), we obtain an estimated value for the photometric B-band 
B/D ratio of 0.19 for ErisBH , which would classify the galaxy 
as a barred Sc spiral. 

We now look at the differences in sfr history between the 
two simulations. The evolution of the sfr of Eris and ErisBH 
if shown in Eigure [T^ (empty and filled squares respectively). 
The symbols of different colors show the sfr averaged over 10 
snapshots (approximately 300 Myr) for the entire box (black) 
and within spheres of different radii from the center. The 
largest differences between the two simulations are in the inner 
few kpc, where the star formation rate can be up to an order 
of magnitude lower in ErisRi/than in Eris. This is reflected in 
the stellar distribution, which is shown in Eigure (8] within few 
kpc the stellar content in ErisBH can be up to 50% lower than 
in Eris, while above 10 kpc there are no striking differences 
between the two simulations. 
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Figure 15. Distribution of “orbital circularity parameter” JcjJz 
for ErisBH (left panel) and Eris (right panel) at 2 ; = 0.3 . The 
distribution including all stars within 20 kpc from the center is given 
by the black lines. The colored lines show the distribution for stars 
within spheres with radial distance from the center as indicated in 
the caption. 


Figure 14. Stellar density profiles of ErisBH and Eris at 2 : = 0 
(solid lines). The two-components Sersic fits are shown by the 
dashed and dotted curves, for the bulge and disk, respectively. The 
values of the best-fit parameters are given in Table 14.21 


Apart for being generally lower in the central regions, the 
sfr in EiisBH has also a different redshift evolution: while in 
Eris the sfr gradually decreased after z = 1, in EmsBH we 
see still some small bursts of star formation at more recent 
times. There is particularly a new increase of star formation 
at z < 0.3, after the appearance of the strong nuclear bar 
visible at z = 0. It is likely that the bar is driving a new 
inflow of gas to the center (see also the gas distribution in the 
central kpc in the left panel of Figure [8} which leads to a new 
episode of star formation. 

The presence of a young generation of stars in the very 
center of ErisBH can be seen also in Figure 1131 where we 
show the age distribution of the stars within different radial 
shells from the center. The very central stars (within 300 pc) 
of EvisBH (solid histogram on top panel) are primarily quite 
old, except for a population of young stars {age < 2 Gyr) 
originating from the new episode of star formation at 2 ; < 0.3 
which takes place during the formation of the bar. In Eris 
(dotted histogram) the age distribution in the inner region is 
different: there is a large population of stars born during the 
last merger at 2 : ~ 1 {age ~ 7 Gyr), and a slowly decreasing 
distribution at younger ages (see the corresponding increase in 
the sfr in the inner region around 2 ; ~ 1 in Figure [T^ followed 
by a gradual decrease of the sfr). 

The role of bars in bringing a population of young 
stars within the bulges o f disk galaxies has been studied by 
(jGoelho Gadot^ l201lh . where the authors speak about a 
“rejuvenation” of bulges due to the presence of bars. 

On top of causing a new inflow of gas which leads to the 
new burst of star formation, the bar has also a dynamical ef¬ 
fect in the distribution of angular momentum in the center of 
the galaxy. Figure [15] shows the distribution of the “orbital cir¬ 
cularity parameter” Jc/ jJB for Eris(left) and Eris (right) 
at 2 ; = 0 for all stars within 20 kpc from the center (black 


curves) and for stars within different radii from the center (col¬ 
ored lines). Starting from the very central kpc (red curves), 
there is a striking difference between EiisBH and Eris: in Eris 
(right panel) there is a narrow peak of stars around Jcj Jz — 0, 
and a significant fraction of higher angular momentum parti¬ 
cles that reach the central region. This is not happening in 
EiisBH , where no high-angular momentum stars are present 
in the center of the galaxy, consistent with the presence of a 
strong bar. The detailed effects of the bar on the dynamical 
properties of the galaxy will be studied in a separate paper. 

Finally, the lower stellar content in the center of EiisBH 
results to a lower rotation curve, with respect to Eris, in the 
inner 10 kpc (Figure (Tb]). A value around 190 — 200 km s“^ 
for the peak of the circular velocity in the inner region of 
the galaxy is consistent with the values estimated for the 
Milky Way by (jPortail et ahllioiBl h who constructed dynam¬ 
ical models of the bulge of the Galaxy using the 3D density 
measurements of Red Glump Giants and kinematic data from 
the BRAVA survey. At large distances, on the other side, the 
rotation curves of the two simu lations are almos t identical, 
and consistent with the data of IXue et al.l (|2008l h who esti¬ 
mated the Milky Way’s rotation curve using the kinematics of 
blue horizontal-branch stars. 

To conclude, despite being quantitatively modest, the 
feedback from the central black hole of EiisBH leads to small 
changes in the gaseous and stellar distribution in the inner 
region of the galaxy, when compared to the Eris simulation 
(where AGN feedback was neglected). Those small variations 
are then the origin of important secondary effects, particularly 
in the stability of the galaxy. EiisBH , in fact, experiences a 
strong instability at 2 : ~ 0.3, that leads to the formation of a 
new large bar. This event causes the depletion of gas in the 
inner few kiloparsecs, and further growth of the bulge is in¬ 
hibited. By 2 ; = 0, the central galaxy of EiisBH has a larger 
disk and a smaller bulge, with a box/peanut morphology, with 
respect to the Eris galaxy. 


5 SUMMARY AND CONCLUSIONS 


® The orbital circularity parameter is used to kinematically sepa¬ 
rate stars belonging to the disk and stars belonging to the bulge. It 
is defined as Jcj Jz = Jz/jcircj where jz is the angular momentum 
of each star in the z-direction (i.e. the direction of J), and jcirc is 
the angular momentum expected for a circular orbit at the same 
radius. 


In this paper we presented EiisBH , a new zoom-in cosmo¬ 
logical simulation where we followed the evolution of a Milky 
Way-size dark matter halo with its baryon content from 2 ; ~ 90 
to the present tim e. EiisBH is a twin simulation of Eris 
(jGuedes et al.l l201lh , with which it shares the same initial 
conditions, simulated physical processes and the resolution of 


© 0000 RAS, MNRAS 000 , 000-000 














14 



Afstar 

Mgas 

fb 

Enax 

Moisk 

MBulge 

B/T 

-Pdisk 

-Pbulge 

Us (Bulge) 

ErisBH 

3.2 X lO^o 

5.9 X 10^° 

0.II8 

193.2 

2. X I0^° 

8.4 X 10^ 

0.29 

2.9 

0.4 

0.9 

Eris 

3.9 X IQio 

5.6 X IQio 

0.I2I 

239.2 

1.9 X IQio 

1.5 X IQio 

0.44 

2.3 

0.3 

1.2 


Table 2. Summary of global properties of KrisBH and Eris at 2 ; = 0. The stellar and gas masses in columns 1 and 2 refer to the total stellar 
and gas masses within the halo, also used to calculate the baryon fraction given in column 3. The disk and bulge masses have been obtained 
fitting a double Sersic profile to the face-on surface density profile, with scale lengths given in columns 8 and 9. The last column gives the 
Sersic index for the bulge, while the Sersic index for the disk has been fixed to unity. Masses are given in solar masses, distances in kpc and 
velocities in km s“^. 


120 pc, but, differently from Eris, it includes also prescriptions 
for the formation and evolution of massive black holes. 

Four black holes are inserted during the simulation run in 
the most massive galaxies that also contain a large gaseous 
high-density region. The largest black hole, located in the 
largest galaxy of the simulated volume, is seeded at z ~ 8.5 
with an initial mass of about 9 x 10^ M©. 

In this first paper on ErisBH we focused on studying the 
properties of this central black hole, its growth, and its relation 
and influence on the host galaxy. Our main conclusions are: 

• During its evolution, very little gas reaches the center 
of the main galaxy, thus little fuel is available for feeding 
the main black hole, whose growth is generally highly sub- 
Eddington. Ending up at 2 ; = 0 with a mass of 2.6 million 
Solar masses, only ~ 3 times its initial value, the black hole 
grows primarily (about a factor of two) through mergers with 
the infalling black holes originally hosted by satellite galaxies; 

• Given the little gas supply, the value of the initial seed 
mass is rather important. Generalizing our results, we expect 
Milky Way-size galaxies to be hosting an intermediate-mass 
black hole already at z ~ 8. This could be possible either via 
a direct collapse in the protogalaxy, or via the uninterrupted 
growth close to the Eddington rate of a Pop III star remnant. 

• The final black hole sits on the Mbu—cf and Mbh — Msuige 
relations in a location close to the one of SgrA*. Gonsistent 
with observational results for pseudobulges, EiisBH is close 
to the extrapolation at lower masses of the Mbh — cr relation 
established by more massive galaxies, while it is almost an 
order of magnitude below the Mbh — Mnuige relation; 

• Given the little growth by accretion of the central black 
hole, the AGN feedback on the galaxy is quite modest. Only 
the very central region of the galaxy (within ~ I — 2 kpc) 
seems to be affected by the presence of the black hole, with 
the consequence that the gas content and the star formation 
rates in the central region of the galaxy are lower than the 
ones found in Eris; 

• Because of its smaller bulge and the lower concentration 
of gas in the center, the disk of EiisBH is more unstable. The 
last instability event, at z ~ 0.3, causes the formation of a bar 
of about 3 kpc in radius and a burst of star formation in the 
very central few hundred parsecs; 

• At z = 0, the disc of EiisBH is slightly larger than the 
one of Eris, while its bulge is almost a factor of two smaller. 
Gonsequently, the B/D ratio of ErisBHis about half of the one 
of Eris and its rotation curve hatter. Moreover, the bulge of 
EvisBHis characterized by the box/peanut morphology typical 
of barred galaxies, and, at its very center, it hosts a population 
of very young stars. 

In summary, as expected, progenitors of late-type galaxies 
seem to not be favourable sites for efficient black hole growth. 



r [kpc] 

Figure 16. Rotation curves for ErisBH (solid lines) and Eris (dot¬ 
ted lines) in the inner 60 kpc. Total vc (black curves), for dark 
matter only (blue curves), for stars only ( orange curves), f or gas 
only (red curves). The data points are from IXne et al.l (|2008h . 


probably due to their rather quiet merger history. While we 
do no witness any dramatic episode of AGN feedback, we hnd, 
however, that the weak, but circumscribed, energy input from 
the black hole helps the global stellar feedback to inhibit the 
growth of the bulge, which has important consequences on the 
global stability of the disk. Modest-size black holes, though, 
might play an important role in the dynamical evolution of 
late-type spirals. 
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